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{S) Abstract 

W) We propose measurement-based quantum computation (MBQC) as a quantum mechanical 

toy model for spacetime. Within this framework, we discuss the constraints on possible temporal 
orders enforced by certain symmetries present in every MBQC. We provide a classification for 
Q\ all MBQC temporal relations compatible with a given initial quantum state and measurement 

CN| setting, in terms of a matroid. Further, we find a symmetry transformation related to local com- 

plementation that leaves the temporal relations invariant. After light cones and closed time- like 
curves have previously been found to have MBQC counterparts, we identify event horizons as a 
Q_i third piece of the phenomenology of General Relativity that has an analogue in MBQC. 

c 

1 Introduction 

i_i Unifying quantum mechanics with general relativity 016] is a major open problem in physics. 

After the unification of quantum mechanics and special relativity has been accomplished 18] , in the 
resulting relativistic quantum field theories [9] - |17j , spacetime and quantum-mechanical degrees of 
freedom are still treated on a separate footing. For example, the electromagnetic, weak and strong 
forces are described as being mediated by particles, namely photons, W and Z-bosons, and gluons. 
An analogous consistent formulation of gravity has so far not succeeded. In order to overcome the 
schism between the theory of gravity on one side and quantum mechanics on the other, one may 
therefore ask: Is there an element in the structure of quantum mechanics that forces spacetime into 
existence? I.e., is spacetime an emergent phenomenon, predicted by quantum mechanics? 
• • In the present paper, we approach this question in a restricted and controlled setting, with 

. ^ emphasis on the notion of time. Namely, we study a scheme of quantum computation [18], so-called 

'measurement-based' or 'one-way' quantum computation (MBQC) [19J as a quantum-mechanical 
toy model for spacetime. In MBQC, the process of computation is driven by local measurements on a 
suitably chosen entangled state. 'Time' in deterministic MBQC is a partial ordering of measurement 
events, which is compatible with the notion of time in the theory of relativity. An external and 
continuously increasing clock time is not required. The possible temporal orders are severely 
constrained by symmetries present in every MBQC. The exploration of these symmetry-enforced 
constraints is a central theme in this paper. 

A fundamental property of quantum mechanics is that measurement outcomes — as opposed to 
measurement bases — cannot be chosen. They are in general random. MBQC, which is driven by 
measurements (no unitary evolution ever takes place in this computational scheme), were entirely 
useless if this randomness could not be prevented from creeping into the logical processing. But it 
can: all that is needed is the adjustment of measurement bases according to measurement outcomes 
obtained at other qubit locations before. One can turn this mechanism for adjusting measurement 
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bases around, and say that it is the inherently probabilistic nature of quantum measurement which 
determines the 'before' in MBQC. Indeed, without even knowing what is being computed in a given 
MBQC, the requirement of preventing the randomness of quantum mechanical measurement from 
affecting the logical processing imposes stringent constraints on the possible temporal orders of 
measurement events. It is because of this link between temporal order and quantum measurement 
that we propose MBQC as a quantum-mechanical toy model for spacetime. 

In the context of the present work, it has previously been shown for the case the initial state 
in MBQC is a stabilizer state [20] known in a suitable local basis, that the entire temporal order 
of measurements is fully specified if the set I of first-measurable and the set O of last-measurable 
qubits are known [2Tj . [22]. But not every pair /, O leads to a possible temporal order. A condition 
on admissible pairs /, O has been given in [23]. On a different route, an MBQC-analogue of closed 
time-like curves (CTCs) has been identified [Mj, via post-selection models [25], [26J for CTCs in 
quantum circuits. Certain solutions of the Einstein equations permit closed time-like curves, and 
if MBQC is to serve as a toy model for a quantum spacetime it is a bonus that it accommodates 
this part of the phenomenology of general relativity. 

The results of the present paper are three- fold. First, we identify a gauge transformation leaving 
the classical output of all MBQCs invariant, and causing additional terms in the classical processing 
relations for measurement outcomes. We show that the resulting generalized processing relations 
fully specify the resource state and measurement setting, up to equivalence. In the reverse direction, 
we provide a classification for all MBQC temporal relations compatible with a given resource 
stabilizer state and set of measurement bases. Specifically, we introduce a matroid capturing all 
information of the resource state and measurement bases relevant for temporal order, and show 
that the classical processing relations — including temporal order — are in one-to-one correspondence 
with the bases of this matroid. Second, we identify a further group of symmetry transformations 
generally changing MBQCs but preserving temporal order. The possible temporal orders in MBQC 
thus arise as representations of this group. These symmetry transformations turn out to be related 
to local complementation [27J-[29j. Third, after light cones and closed time-like curves, we identify 
event horizons as a third piece of GR phenomenology that has a counterpart in MBQC. 

2 Background 

In this section we review some basic facts about measurement-based quantum computation, es- 
sential definitions for the discussion of MBQC temporal order, as well as previous work in this 
area. The scheme of MBQC itself, with a proof of its computational universality, is not reviewed 
here since various articles on this subject exist in the literature [T5], [30] - [34J Also, we require 
familiarity with the stabilizer formalism |20| . 

2.1 MBQC and cluster states 

In MBQC, the process of computation is driven by local (=l-qubit) measurements on an initial 
highly entangled state, generally taken to be a so-called cluster state. The local measurements can 
only reduce entanglement, and therefore all entanglement needed for the computation must come 
from the initial state. For this reason, the initial state is often called the 'resource stateQ 

The computational power of a given MBQC strongly depends on the choice of the initial resource 
state. For example, a local resource state has obviously no computational power. Other states 
may be used for a restricted class of computations. Two-or-higher dimensional cluster states of 
unbounded size have the property that they enable universal quantum computation. That is, any 

1 Large entanglement is a necessary [351 136] but not sufficient condition for the usefulness of a quantum state as 
resource in MBQC. Somewhat paradoxically, quantum states exist that are too entangled to be useful [37] ■ 
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quantum computation can be realized on such a state, by suitable choice of the local measurement 
bases. 

We now define graph states, to be used later on, and cluster states as a subclass thereof. 

Definition 1 (Graph states and cluster states.). Be G a graph with vertex set V(G) and edge set 
E(G), such that there is one qubit for each vertex a G V(G). Then, the graph state \G) is the 
unique (up to global phase) joint eigenstate, \G) = K a \ G) , for all a G V(G), of the operators 

K a = o^ (g) ai b \ (1) 
b\(a,b)eE(G) 

The cluster state \<&c) * s a graph state where the corresponding graph is a d- dimensional lattice C 
In the standard scheme |19j . the measured observables are of the form 

O a [q a ] = coB<pa4 tt) + ("I) 9 * sin ^4°), (2) 

with a the qubit to which the measurement is applied, and q a G 7L<i depending on outcomes of 
(earlier) measurements on other qubit locations. 

2.2 Temporal order in MBQC 

As noted above, the temporal order of measurement in MBQC is a consequence of the randomness 
of measurement outcomes. By adjusting measurement bases according to measurement outcomes 
obtained on other qubits, this randomness inherent in quantum mechanical measurement can be 
kept from creeping into the logical processing [19J. If the measurement outcome of qubit a influences 
the choice of measurement basis for qubit b, clearly, qubit a must be measured before qubit b can. 
This is how a temporal order among the measurement events arises. 

We now generalize the above scenario of measuring observables of form Eq. ^ on cluster 
states. Namely, we now consider general stabilizer states \^f) as resources, with supp(|\I/)) = Q 
and stabilizer group Furthermore, we generalize the local observables whose measurement 

drives the computation from Eq. ^ to 

O a [q a ) = cos ip a af + (-l)*'sin ¥ > a <7$, Va G $7, (3) 

with 00, <7 S ^ G {X, Y, Z}. Therein, the measurement angles <p a are in the range — ir/2 < (p a < 
7r/2, and q a G Z2 may depend on measurement outcomes from several (other) qubits in Q. 

Definition 2 (Measurement plane). For every qubit a G supp(\^/)), the measurement plane at a is 
the ordered pair \ai ° , o~g)\. 

We define a third Pauli operator, a s = ia S( ^a^. As we will see shortly, the Pauli operators a a, 
and a s are useful because of the relations 

a^O[q]al = 0[qel], a s O{q]at = -0[q\. (4) 

The basic mechanism of accounting for an "undesired" measurement outcome is the following. Sup- 
pose on some qubit a G f2, instead of the "desired" post-measurement state \<p a } a the "undesired" 
post-measurement state \<Pa)a has been obtained. The goal is to get the computation back on 
track by only adjusting the subsequent measurements. To do that, we require a stabilizer operator 
K(a) G 5(1^)) with the following properties |21j: (1) K(a) has support only on a and the yet un- 
measured qubits, and (2) K(a)\ a = ai a \ Recall that o s \tp a ) = \<Pa) fc> r the eigenstates \<£ a (qa, s a )) 
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of the local measured observable 0[q a ], c.f. Eq. Q. Denote by V(a) and J-"(a) the past and future 
of a, respectively. Then, 

{v(a)(<floc\ ® a(<fa\) I*) = {v(a)(floc\ ® a(^|) f (o)|*) ^ 

= (p(a)(<Aoc| ® 0(^0 1) ^(«) |.F(a) 1^) 

Therein, the first equality follows from K{a) G <S(|\I/)), and the second from the above properties 
(1) and (2). 

Since the overlaps between local states (representing the local measurements) with the resource 
state \^f) contain all information about the computation, we thus find that we can correct for 
"undesired" outcomes by (a) adjusting measurement bases of future measurements (caused by 
tensor factors as in K{a) \j^t a ) ) an d (b) re-interpretation of measurement outcomes (caused by 
tensor factors a s in K[a) \j^( a ))- 

Example: Consider MBQC on a cluster state of three qubits on a line, each measured in 
the [a x , cjjj-plane. That is, for all three qubits a s = Z and as = X. (Here and from now on, we 
use the shorthand X = a x , Y = a y and Z = a z .) The stabilizer generators of are 

K x = X 1 ®Z 2 ®I 3 = ®ai 2) ®/( 3 ), 

K 2 = Z l ®X 2 ®Z 3 = al 1] ®a^ ] ® af\ (6) 
K 3 = h®Z 2 ®X 3 = /« <8 a{ 2) <8 af . 

When the three cluster qubits are measured in the order 1 -< 2 -< 3, the corresponding quantum 
circuit is EH 



+>-U z — U x — U z -\\X (7) 



Using Eq. ([5|, here we show that if qubits 1, 2, 3 are measured in the order 1 -< 2 -< 3, then the 
randomness of the measurement outcomes on qubits 1 and 2 can be corrected for. First we consider 
the stabilizer operator K 2 = K(l). If inserted in the state overlap of Eq. (pfo, 

the measurement outcome of qubits 1 and 3 is flipped, as well as the measurement basis at quBit 
2. Since qubits 2 and 3 are yet unmeasured when qubit 1 is measured, this is a valid correction 
operation for qubit 1; hence the notation K(l). Similarly, K 3 = 1^ (8) a^ ® of 1 =: K(2) can 
be used as correction operation for qubit 2. K(2) flips the measurement outcome of qubit 2 and 
the measurement basis of qubit 3. Since qubit 3 is yet unmeasured when qubit 2 is measured, this 
corresponds to a valid correction operation. 

The above argument also works in reverse. If the correction operations K(l) and K(2) are 
used, then 1 -< 2 -< 3 follows. K{1) implies that the measurement basis of qubit 2 depends on the 
measurement outcome at qubit 1, hence 1 -< 2. Note that 1 -< 3 does not yet follow! K(l) does 
not affect the measurement basis at qubit 3. Only the meaning of the eigenstates is interchanged, 
which by itself does not require qubit 3 to be measured after qubit 1. The interpretation of the 
measurement outcome may take place long after the measurement itself has taken place. 

K 3 = K(2) implies that the measurement basis of qubit 3 depends on the measurement outcome 
of qubit 2, and hence 2 -< 3. Both relations taken together yield 1 -< 2 -< 3. 

From the equivalence with the circuit of Eq. Q one would expect one bit of classical output. 
Indeed, if no correction operations need to be used, the eigenvalue measured at the output of the 
circuit corresponds to the eigenvalue A3 measured on qubit 3 of the cluster. Now recall that K(l), 
applied conditioned upon Ai = —1, flips A3. Therefore, with or without corrections, the eigenvalue 
measured in the circuit Eq. ([7]) equals A1A3. Or, in binary notation Ai = (— l) Sl , A3 = ( — 1) S3 , the 
single bit of classical output takes the value s\ + S3 mod 2. 

Note that we have not made any use of K\ = a} (8) in the above argument. Still, K\ has 
a role to play, as we discuss in Section |4| 
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2.3 Influence matrix, forward and backward cones 



To counteract the randomness of measurement outcomes, two measurement settings (i.e., bases) 
per qubit suffice, which may be labeled by q a = and q a = 1, respectively, for each qubit a. The 
measurement settings may collectively be described by a binary vector q, [q] ffl = q a for all a£(l, 
and the measurement outcomes by a binary vector s, [s] a = s a , for all a G O. It turns out that the 
relation between measurement bases q and measurement outcomes s is linear [38J, 



q = Ts mod 2, 



(8) 



with T a binary matrix. We call T the influence matrix. 

The set of all qubits b whose measurement basis must be adjusted according to the measurement 
outcome on a is denoted as the forward cone of a qubit a. Similarly, the backward cone of a qubit 
b is the set of all those qubits a whose measurement outcome influence the measurement basis at 
b. More formally, with Eq. Q, 

Definition 3 (Forward and backward cones). For any a G 0, the forward cone fc(a) is given by 

fc(a) := {b G fi| dq b /ds a = 1} . (9) 

For any b G f2 the backward cone bc(b) is given by 

bc{b) := {a G n\ dq b /ds a = 1} . (10) 

We denote the characteristic vectors of fc(a) and bc(b) by fc(a) and bc(6), respectively. Then, 
the influence matrix T takes the form 



/ 



bc(l) ) \ 
bc(2) ) 



V ( bc(n) ) J 



( \ 



fcfr 



/ \ 

fc(2) 



VV / V / 



fc(n 

V J J 



(11) 



The influence matrix T generates a temporal relation among the measurement events under tran- 
sitivity. We say a <b (a precedes b) if 6 G fc{a). A priori, it is not forbidden that for two qubits 
a,b G O,, a -< b and b -< a. However, such a computation could not be run deterministically in 
a world like ours where time progresses linearly. An MBQC is deterministically runnable if the 
temporal relation "-<" between the measurement events is a strict partial order. 

Definition 4 (Strict partial order). A strict partial order is a relation among the elements a G £1 
with the following properties 



a -ft a, 

a -< b ==> b -fi a, 

a -< b,b ~< c ==> a ~< c, 



Va G Q, 
Va, b G 0, 
Va, 6, c G £1. 



(irreflexivity) 
(antisymmetry) 
(transitivity) 



(12) 



Definition 5 (Input and output sets.). For a given MBQC, the input set I CS1 is the set of qubits 
whose backward cones are empty, I = {a G fi| bc(a) = 0}. The output set O Q £1 is the set of qubits 
whose forward cones are empty, O = {b G O, fc{b) = 0}. 

That is, with respect to a given temporal relation among the measurement events, / is the 
maximal set of qubits which can be measured first, and O is the maximal set of qubits which can 
be measured last. 
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Regarding the computational output, for the purpose of this paper we are exclusively interested 
in MBQCs for which the computational result is a classical bit string. That is, every qubit in 
is measured. Then, as a consequence of the randomness of individual measurement outcomes, the 
classical output o of an MBQC is given by correlations among measurement outcomes. Again, the 
relation between classical output and measurement outcomes is linear, 

o = Zs mod 2, (13) 

for a suitable binary matrix Z. 

Example: To illustrate the above notions, we briefly return to the three-qubit cluster state 



example of Section |2.2[ From the previous discussion we find that 

mod 2. (14) 





Therefore, fc(l) = {2}, /c(2) = {3}, /c(3) = and 6c(l) = 0, 6c(2) = {1}, 6c(3) = {2}. Hence, 
I = {1} and O = {3}. Also, 1^2 and 2 -< 3. The latter two relations generate a third under 
transitivity, namely 1 -< 3. Asymmetry and irreflexivity are obeyed in this example. Regarding the 



single bit of output in this computation, the matrix Z of Eq. (13) is Z = (1 01 



2.4 Brief review of prior work on MBQC temporal order 

MBQC temporal order as an (almost) emergent phenomenon. In [21] the following question is asked: 
"Given graph G and the set E of measurement planes for all vertices, can the temporal order of 
measurements in MBQC with a graph state \G) be uniquely reconstructed from this information?" 
The graph G and the measurement planes E are undirected objects. Thus, if the answer to this 
question was yes, then temporal order in MBQC were truly emergent. 

However, it turns out that the pair G, E does not specify the temporal order of measurements 
in MBQC uniquely; there are in general a number of consistent temporal orders respecting the 
requirement that the randomness of measurement outcomes should not affect the logical processing. 
One may then ask how constraining on temporal order this requirement actually is. To this question, 
the following answer is provided by [21] : If in addition to G and E the set I of first-measurable and 
the set O of last-measurable qubits is known, then the complete temporal order (if existing) can 
be uniquely reconstructed from this information. Thus, MBQC temporal order is not emergent in 
the strict sense; a seed /, O must be provided in addition to G and E, and the complete temporal 
order then follows. 

But not every pair /, O will lead to a consistent temporal order. The question that now arises 
is which pairs /, O do. For the case where the stabilizer resource state is a graph state and all 
qubits are measured in the [X, y]-plan^] then the answer to this question is given in [23j. Denote 
by Aq the adjacency matrix of the graph G describing the resource state |G), and by Ag\o c xI c 
the submatrix of Aq where the rows are restricted to O c := D\0 and the columns are restricted to 
I c . Then, the pair /, O leads to a partial order of measurement events in MBQC iff there exists a 
matrix T such that Ag\o c xI c T = I, and T is free of cycles (that is T aa = 0, Va, T^T^ = 0,Va, b, 
T a bTb c T ca = 0, Va, b, c, etc). The resulting temporal order is generated by T under transitivity. 

MBQC and closed time-like curves. In [21] it is shown that MBQC encompasses the post-selection 
model of closed time-like curves (CTC's) proposed by Bennett, Schumacher |25| and Svetlichny 
|26| . The CTCs arise from circuits such as the one displayed in Fig. [T^,, translated into MBQC. 

2 The former condition alone is not a restriction, since all stabilizer states are local Clifford equivalent to graph 
states [39], but both conditions jointly are. 
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a) b) 




Figure 1: Closed time-like curves in MBQC. a) Bennett, Schumacher and Svetlichny's post-selection 
model [251 [26] of CTCs (left: circuit with wires 'going backwards in time', right: implementation 
thereof using teleportation and post-selection), b) Nested forward cones in the MBQC equivalent 
of the teleportation circuit in (a). 



The result are forward cones with the property b G fc(a) Ao £ fc(b), for two suitably chosen qubits 
a, b G S7; see Fig. [Tja. Such nested forward cones are an obstruction to deterministic runnability of 
MBQC, but mimic closed time-like curves of General Relativity in the MBQC setting. 



3 Outline 

The remainder of the present paper is organized as follows. In Section|4]we show that every quantum 
computation — both in the circuit and the measurement-based model — transforms under a gauge 
transformation which affects the sequence of gates or measurements, respectively, but keeps the 
probability distribution of computational outputs invariant. For MBQC, this gauge transformation 
has the effect of causing extra terms in the classical processing relations Eq. Q, (13). 



In Sections [5] and [6] the consequences of these gauge transformations for MBQC are explored. 
It is already known [21J that the classical processing relations (containing the temporal order 
of measurement events) are not independent of the resource state and the set of measurement 
planes. For the extended processing relations — which now contain more information — this mutual 
dependence becomes stronger. In fact, the classical processing relations determine the resource 
state and set of measurement planes, c.f. Theorem [4] in Section [5j 

In Section [7] we introduce a second symmetry transformation acting on MBQCs which we call 
'flipping of the measurement plane'; see Fig (2ja. As introduced in Section [2j the measurement 
plane [a^,a S( j,] is an ordered pair of Pauli operators. Or is it? What happens under the exchange 
<70 { — > <r s 0? (Under this transformation, the set of measurable observables in the plane is mapped 
back onto itself, hence the name for the transformation.) Denote by a the qubit whose measurement 
plane is flipped. It turns out that whenever the measurement basis at a does not depend on the 
measurement outcome at a, i.e., in the absence of a self-loop at a, flipping of the measurement 
plane at a induces a transformation on the influence matrix T which leaves the temporal order 
generated by T invariant. The effect of flipping measurement planes on the influence matrix T is 
closely related to the operation of local complementation |27]-|29j on graphs. 

In the present context, we consider MBQC as a toy model for spacetime and therefore admit 
temporal relations with closed time-like curves. However, closed time-like curves are an obstacle to 
deterministic MBQC. In Section [8] we show that every MBQC solving a problem in the complexity 
class NP can, without reducing its probability of success, be transformed into a deterministically 
runnable MBQC, i.e., into an MBQC whose temporal relation of measurement events is a strict 
partial order. This transformation does not guard against the inefficiencies of post-selection, but 
restores deterministic runnability. 

In Section [9] we conclude and point out open questions. 
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a) b) 




Figure 2: Two symmetry transformations, a) Gauge transformation in the circuit model. For any 
logical qubit I, an identity I = &x &x is inserted into the circuit next to the input. The left a x is 
propagated backwards in time, and absorbed by the input state |+). The right o x is propagated 
forward in time, flipping rotation angles and, potentially, measurement outcomes in the passing, b) 
Flipping a measurement plane in MBQC. In the measurement plane [cr^,, cx,,^], the Pauli operator 
(70 is distinguished over <r S( ^ because the rule for adjusting a local observable O a for measurement 
is O a [q a = 1] = <J^O a {q a = 0](f^ ) ■ If T aa = (with T the influence matrix) then the flip 
Cj*' < — > is a symmetry transformation for the given MBQC. 



4 Gauge degrees of freedom 

Here we introduce the notion of "gauge transformations" acting on a given quantum computation. 
These transformations exist for both the circuit model and MBQC. 



4.1 Gauge transformations in the circuit model 

To obtain an intuition for the gauge transformations introduced here, it is instructive to first inspect 
them in the circuit model. Specifically, we consider a quantum circuit which consists of (1) the 
preparation of the quantum register in the initial state ®" =1 (2) unitary evolution composed 
of, say, CNOT gates and one-qubit rotations about the X- and Z-axes, and (3) local measurements 
for readout. Such a circuit is displayed in Fig. [2] above. Then, into every qubit line individually, we 
may insert an identity / = cr x a x next to the input; See Fig. [2] The left a x is propagated backwards 
in time until absorbed by the input state |+). The right a x is propagated forward in time, flipping 
rotation angles and readout-measurement outcomes in the passing. 

This transformation is an equivalence transformation, since it is caused by the insertion of an 
identity gate into the circuit of Fig. [2j The reason we call it a gauge transformation is that it 
changes reference frames for certain rotation angles. Specifically, for the rotation gates next to 
each input qubit we can individually choose our convention for which rotation angles are called 
positive or negative, respectively. Once those signs are fixed on the input side, they are fixed 
throughout the circuit. 



4.2 The gauge transformations in the measurement-based model 

Translating the above discussed gauge transformations from the circuit model into MBQC it is 
easily seen that the above relations Eq. ([8| and (13) are incomplete. We find the more general 
relations 

q = Ts + Hg mod 2, (15a) 
o = Zs + Rg mod 2, (15b) 



S 



with g a choice of gauge. 

Under a change of g, the measurement bases in a particular MBQC change, but the probability 
distribution for the classical output values remains unchanged. Of course, knowledge of the g-de- 
pendent extra parts in Eq. ( 15a|15b ) is not necessary to run any given MBQC, since g = is always 
a valid choice. However, the presence of the extra terms in the processing relations strengthens 
their interdependence with the resource state, which is the reason why we discuss them here. 

We may now want to derive the generalized processing relations Eq. ( 15a|15b ) directly in MBQC, 
without reference to the circuit model. Before discussing the general case, we return to the specific 



example of the three-qubit cluster state in Section 2.2 

Example: Consider the product K\K^ = of K\, K3 in Eq. (6). When used in Eq. (5), 

the effect of this stabilizer element is to flip the measurement bases of qubits 1 and 3. Hence, the 



relation Eq. ( 14 ) generalizes to 















Q2 


H 


1 








Q3 J 




\o 


1 







gi mod 2. 



(16) 



This is precisely what one would expect from the insertion of a x a x next to the input of the equivalent 
quantum circuit in Eq. ([7]). 

We now turn to the general case. Via Eq. ([5]), the stabilizer group <S(|\1/)) acts on s and q. 
Denote the post-measurement state of qubit a G CI by \s a ,q a ) a , with s a the measurement outcome 
and q a specifying the measurement basis. Then, as in Eq. ([5]), 



, a (s a ,q a \ I I*) = I QS)a(s a ,qa\ I K\V) 



\ a (s a ,q a \K I*), MKzSm). (17) 



Since a s \s,q) = |s® 1, q) and a,p\s, q) = |s, q® 1), for a stabilizer element K = (S) a en( cr s ) Va {& t 
the action of Gk on s, q is 



G K ■ 



s — > s + v mod 2, 
q — > q + w mod 2. 



(18) 



Again, nothing changes by the insertion of a stabilizer (identity) operator into the state overlap of 
Eq. ([5]), and transformations Gk are therefore equivalence transformations. They can be used to 
constrain the possible temporal orders in MBQC, as we discuss explicitly in Section [6j 



4.3 Correction and gauge operations in the stabilizer formalism 

To state and prove our results on MBQC temporal order, we need to make a few more definitions. 
First, let us emphasize that in this paper we allow closed time like curves (CTC) in the computation; 
we are interested in classifying all transitive temporal relations consistent with a given resource 
stabilizer state and set of measurement planes. Therefore, we do not restrict to partial orders per 
se, and drop the conditions of irreflexivity and antisymmetry on the temporal relations. 

It turns out that the possible temporal relations, and indeed the classical processing relations 



Eq. (15a), (15b), can be parametrized by two subsets of Cl, namely the computational output set 
Ocomp and the gauge input set /gauge- We define these sets next. 

We say that the measurement outcome s a of qubit a is corrected in a given MBQC, if by insertion 



of a suitable stabilizer operator K in Eq. (17) equivalence of s a = 1 with the reference outcome 
s a = is established, at the cost of adjustment of measurement bases and/or re-interpretation of 
measurement outcomes on other qubits. 
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We observe that if for a given MBQC all measurement outcomes s a , o£fi, can be corrected 
then these measurement outcomes contain no information and no linear combination of them is 
worth outputting. Thus, in general there will be a set of qubits whose measurement outcomes are 
not corrected. 

Definition 6 (Computational output set). For a given MBQC, the computational output set 
Ocomp C O is the set of qubits whose measurement outcomes are not corrected. 

The correction operations for the qubits in a £ (O com p) c 5 c.f. Eq. ([5]), are each implemented 
by correction operators K(a) £ For each a £ (O CO mp) c , K(a) has the property that it flips 

the measurement outcome s a at qubit a, but does not flip the measurement outcome on any other 
qubit in (0 C omp) c - I n this way, it is ensured that the correction operation is for qubit a individually. 

Definition 7 (Correction operator). For an MBQC with a given stabilizer resource state \^f), 
fixed set S of measurement planes and computational output set O comp , for each a £ (O comp ) c the 
corresponding correction operator K(a) G 5(1^)) is a Pauli operator satisfying the conditions 



K(a)\ a G {ai a \aip, 

K(a)\ b £ {/ ( Vi ft) }, V6 G [QcompYXa. 



(19) 



The correction operator K(a) can be used to correct an "undesired" measurement outcome at 
qubit a, c.f. Eq. (5). If, for a given operator K (a), K(a)\{, £ {a, , , cr5} then, by Eq. (5), the 
at qubit b depends on the measurement outcome for qubit a. In terms of the 



measurement basis 
influence matrix T. 



K(a)\ b £{a { '\afh^[T] ba = l, Va £ (O comp ) c , b £ Q. (20) 



Note that in the first line of Eq. (19) we allow K(a)\ a = a^f) only because we are admitting closed 



time-like curves in the present discussion. K(a)\ a = a^l means that the measurement basis for 
qubit a depends on the outcome s a of the measurement of qubit a. This amounts to a closed 
time-like curve only involving qubit a (self-loop) and is an obstacle to deterministic runnability. 
Because the qubits a £ O com p have no correction operations, T ba = for all b £ £1, and thus 

Ocomp C O. (21) 

Thus, C omp and the correction operators {K(a),a £ (0 C omp) c } completely determine the influence 
matrix T and hence the temporal relation among the measurements. 



Conversely, if O CO mp and {K(a),a £ (O CO mp) c } are unknown, then the constraints Eq. (19) pose 
self-consistency conditions on them. We discuss these conditions further below. 

As we have seen in the concrete three-qubit example above, for a given MBQC the relation 
between the choice q of measurement bases and the measurement outcomes s in general allows for 



an offset term Hg, c.f. Eq. (15a). Thus, the measurement bases for a certain set of qubits can be 



freely chosen until the initially arbitrary g becomes fixed. This observation leads to 

Definition 8 (Gauge input set). For a given MBQC with input set I, the gauge input set I gauge Q I 
is a set of qubits such that for each i £ Igauge the parameter qi specifying the locally measured 
observable Oj[%] can be freely chosen. 

Analogously to the correction operations, there will be gauge operators which implement the 
'corrections' of measurement bases for the qubits in /gauge- The definition of the gauge operators 
ensures that for all i £ /gauge the corresponding Oj can be changed individually without changing 
the others. 
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Definition 9 (Gauge operators). For an MBQC with a given stabilizer resource state \^f), fixed 
set £ of measurement planes, computational output set O comp and gauge input set I ga uge, for each 
i G Igauge, the corresponding gauge operator K{i) G 5(1^)) is a Pauli operator satisfying the 



conditions 



K(i)\i 
K{i)\, 
K(i)\ k 
K(i)\i 



G 
G 



a {i \ 



VfcG((J, 



G {Igauge H {Ocomp) )\^> 



gauge) H {Ocomp) )\^> 



(22) 



VZ G {Igauge H Ocomp 



Like previously for the correction operations, Eq. (22) poses self-consistency condition on the 
possible sets igauge and {if(i),i G /gauge}- 



Lemma 1. T/ie correction operators K{a) of Eq. (19), a G {O com p) c , and the gauge operators 
K{i) of Eq. (22), i G Igauge, are independent. That is, for all sets J C Igauge, L C {O comp ) c with 



J jL V L ^ 0, K(J,L) := n ie j^« n a£ L^(«) ^ ^ (n) - 

Proof of Lemma^ (indirect) Assume that there exists a pair J, L, with J / V L / 0, such 
that if(J,L) = / {f7) . Then, for all 6 G (O comp ) c , K{J,L)\ b = 1. Now^with Eq. (l9], if(6)| 6 G 



{a^exfj}, and 2f(c)| 6 G MV< b >} for all c G (O comp ) c \6. With Eq. fal K{i)\ b G {<ri 6 \ J< 6 >} for 



J>>) 



'sip 



all i G /gauge- Therefore, no other K{-), K{- 



can cancel a ai 



Hence, b G" L for all 6 G (O, 



comp y 



and thus L 



contribution from K{b) to K{J,L). 
By an analogous argument, no K{-) can 



cancel the a, -contribution of K{i) to K{J,L), hence i G" J for all i G /gauge; and J 



Thus, 



# ( J, L) = I (n) => J, L 



Contradiction. Hence, the if (a), if(i) are independent. □ 



Remark: For any given MBQC with fixed temporal relation, the computational output set 



Ocomp is a subset of the output set O, see Eq. (21). But O comp and O are not necessarily equal. 



Likewise, /gauge ^ / by definition, but /gauge and / may not be equal. To illustrate this point, we 
consider the following two examples. 

Example 1. Consider the 3-qubit cluster state of Section 2.2, with measurement planes [-Xyy] 
for all three qubits. We consider the correction operators K{\) = K2, K{2) = K3 for qubits 1 
and 2, and gauge operator K{1) = K1K3. C omp = {3} and /gauge = {1} are then permitted by 



Eqs. (19) and (22), respectively. As was discussed previously for the above choice of correction 
operations, / = {1} and O = {3}. Thus, in the present example /gauge = / and O c 



'comp 



o. 



Example 2. Consider a Greenberger-Horne-Zeilinger state \GHZ) = ( 1 000) + |lll))/\/2 as 
resource state, with all three qubits measured in a basis in the [X, Y] -plane. We use the stabilizer 
elements K{1) := Z\Z% = aids and K{2) := Z2Z3 = ai^a^ as correction operations for qubits 
1 and 2, and K{\) := X1X2X3 = cr^ as gauge operator. Then, the choice / gaU ge = {1} and 

o 



comp 

of a temporarily flat MQC, T 



{3} is admitted by Eqs. (19) and (22), respectively. On the other hand, this is an example 

{1, 2, 3}, and C omp 7^ O, /gauge 7^ /• 



0. Therefore, 1 = 



Lemma 2. For any MBQC on a stabilizer state, \I ga uge\ 21 \^comp\- 



< IO, 



We prove Lemma [2] in Section 4.4 



Definition 10. A pair I gaU ge, O comp is called extremal iff\I, 



1 gauge \ 



IO 



comp I • 



As will become clear in the next section, extremal pairs /gauge; O, 
general pairs, and are not very restrictive (c.f. Theorem [2|. 



comp 



are easier to handle than 



We still need to relate the classical output vector o appearing in Eq. (15b) to the set O comp . 
To this end, we make the following 
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Definition 11 (Optimal classical output). A classical output vector o, with processing relations 
o = Zs + Rg, is optimal iff the following conditions hold 



1. Maximality: Upon left-multiplication by an invertible matrix, Z can be brought into a unique 
normal form 

Z ~ (Z|I) , (23) 
where the column split is between (O com p) c \O cornp , and 



2. Determinism: For the matrix Z in Eq. (23), 



[Z]« = l<=>Jr(;)lie{a« > a®}. (24) 



In other words Eq. (24) informs us that the jth column of Z is simply the restriction of the 

Support Of K(j) to C omp- 

The reason for defining an 'optimal classical output' besides a 'classical output' is the following: 
One could, in principle, run an MBQC perfectly deterministically and then choose o such that 
nothing is outputted at all, or all outputted bits, independent of the measurement angles chosen, 
are zero guaranteed or perfectly random guaranteed. The above definition of an optimal classical 
output eliminates such choices. Maximality says that there is one bit of optimal output per qubit of 
Ocomp- The determinism condition can be understood from the correction procedure explained in 
Eq. ([5]). For a qubit j G (0 C omp) C ! we account for the undesired outcome Sj = 1 by inserting K(j) 



into the state overlap ((piod^) = (<pi oc \K(j)\^). Now consider a qubit i G Ocomp- By Eq. (23), s 
contributes to the output bit Oj, the i-th bit of o. If K(j)\i G {<r« , c^}, then the insertion of K(j) 
into the overlap flips s$, i.e., Si — > Si © 1. This needs to be taken into account when reading out 
Sj. The linear combination Sj © Sj remains unaffected by the correction for Sj. 

We have now made the definitions needed to state and prove our results on temporal order in 



MBQC. In Section 4.4 we establish a normal form for the stabilizer generator matrix of the resource 
state |\&). This normal form is the basis for our results on MBQC temporal order, which are stated 
in Sections [H and HI 

4.4 Normal form of the resource state stabilizer 

Let us briefly review which pieces of information specify an MBQC. A priori, there are four: 
the set of measurement angles, the set S of measurement planes, the resource state \^f) and the 
classical processing relations Eq. ( 15a|15b ). The measurement angles entirely drop out of all our 



considerations about temporal order. Next, we observe that when specifying the measurement 
planes and the resource state separately, we really specify too much. Starting from a given pair 
S of resource state and set of measurement planes, for any local Clifford unitary U, the pair 
U\^),U (£) obtained by applying U to both the measurement planes S and the stabilizer state \^f) 
is again a valid pair, i.e., it consists of a set of a stabilizer state and a set of measurement planes. 
Furthermore, it amounts to exactly the same computation as the original pair. The pair is 
thus redundant. To remove this redundancy, we combine the measurement planes and the stabilizer 
state l^) into the stabilizer generator matrix in the cr^/ovstabilizer basis, 

0(|¥» = (*||5). (25) 

Therein, the columns to the left (right) form the (o~ s -) part of the stabilizer generator matrix. 
<5(| , I / )) comprises all information from the resource state |^) and the set of measurement planes 
£ relevant for the discussion of MBQC. We do not need to know the state and the measurement 
planes separately. 
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Now note that the correction operators K(a), a G (O. 



comp) c , and the gauge operators K(i), 
i G /gauge, are all elements of the stabilizer <S(|\]/)), and, by Lemma[TJ are independent. Thus, they 
either form or can be completed to a set of generators for This observation leads us to the 

following 



Lemma 3. For any MBQC on a stabilizer state with extremal I, 
matrix Q ofS(\^)) can be written in the normal form 



gauge; O com p, 



the generator 



a. 





I gauge 


{.I gauge) {Ocomp 


1 ^ comp 




( 






\ 







T T 


I 


Z T 




i 


H T 









\ 






J 



(26) 



The matrices H, R, T, Z are related to the matrices H, R, T , Z governing the classical processing 
of measurement outcomes in MBQC, q = Ts + Hg mod 2 and o = Zs + Rg mod 2, via 



T 








T 






H 



Z 



Z|J) , R = R. 



(27) 



Proof. By Definition [7j a correction operator K(a) exists for every a G (O c 



'com P ) c . By Eq. @ these 

operators have no support in / D -^gauge and are |(0 C omp) c | in number and must take the following 
form: ( | A | / | B ), where the part is split as /gauge; ^gauge while the cr s -part is split as 
Ocomp, (O comp ) c . The gauge operators are, by definition, of the form: (/C O D), where the 
column splits are as for the correction operators. Since there are |0| generators for the stabilizer 
Q(\^)), of which |(O com p) c | are already accounted for, there can be at most |O comp | independent 
gauge operators, i.e., \L 



augcl 



< o, 



comp | 



/gaugo Ocomp is extremal by assumption, thus these two sets of generators exhaust the stabilizer 



which proves Lemma [2j In the present setting, the pair 
i, thu 

generators and we can write the stabilizer as 



/gauge (/gauge) (O com p) O CO mp 





( 


A 


I 


B 


5(1*}) = 


I 1 


C 





D 



\ 



(28) 



for suitable matrices A, B, C and D. We now need to identify these matrices. By definition, 



'gauge 



C I. Measurement outcomes on qubits a G O, 



comp are not corrected, hence fc(a) = for all 

Then, the influence 



a G O comp , and O comp C O follows from the definition of the output set O 
matrix T takes the form 








T 






(29) 



where the column split is (O, 



comp ) 



C omp and the row split is / ga uge|(/j 



Now consider the 

correction operator K(a) for a G (O comp ) c in the upper part of Q{\^>)) in Eq. (28). We already 



i gauge; 



know from Eq. (20), that the 



part of K{a) is the forward cone of a. Therefore we must have 
Further comparing, with Eq. (24) which states that the restriction of the a s part of the 



A = T 

correction operator K(a) to O comp is the ath column of Z. But this is precisely the ath row of B 
thus we infer that B = Z T . 
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Next, consider row i of the lower part of in Eq. (28k. Row % is (0, .., 0, 1, 0, .., 0|c T | |0|d^). 

The corresponding stabilizer operator K(i), when inserted into the overlap (ip\ oc \^} as in Eq. (|5j), 
flips the measurement basis at qubit i and of qubits I 6 (/gauge ) c with [c]j = 1. It further flips the 
measurement outcomes at qubits m £ O CO mp with [d] m = 1. Therefore, 

H = (—) , with H = C T , and 



x H / 

R = R, with R = D 



We thus arrive at the normal form Eq. (26). □ 



5 Interdependence of resource state and temporal order 



Let us return to our discussion from the beginning of Section |4.4[ on which pieces of information 
are needed to describe an MBQC. At this stage, apart from the set of measurement angles which 
do not enter our discussion, we remain with two pieces of data specifying a given MBQC, namely 



G{\^)) and the processing relations Eq. (15a), (15b). But it doesn't stop there. As the results of 



|21j . |22j . [23] show, the resource state \^), the measurement planes £ and the influence matrix T- 



being part of the classical processing relations Eq. (15a), (15b) — are not independent. Specifically, 
given the sets I of first-measurable and O of last-measurable qubits in addition to and S, the 
temporal order (generated by T) can be worked out completely. 

Here we prove a statement about the interdependence of Q(\^)) and the MBQC classical pro- 
cessing relations which goes the opposite direction. Namely we show that the classical processing 



relations Eq. (15a), (15b) uniquely specify the stabilizer generator matrix Q(\^}), i.e., the pair \^>), 
£ up to equivalence; See Theorem [4] below. Thus, only two pieces of data are needed to specify an 
MBQC that satisfies the determinism constraints, namely the measurement angles and the classical 
processing relations for the measurement outcomes. 

A further question is whether the temporal relations compatible with a resource state |\&) and 
set of measurement planes £ fit into a common framework. In this regard, we show that the classical 
processing relations (containing the temporal order) for MBQC with a fixed resource state and set 
of measurement planes, for extremal pairs / ga uge> O comp , are in one-to-one correspondence with the 
bases of G(\V)), c.f. Theorem [3} 

5.1 Results 

We now present four theorems on the mutual dependence of the resource state and the classical 
processing relations. 

Theorem 1. Consider an MBQC on a stabilizer state \^f), with fixed measurement planes and an 
extremal pair of gauge input set I gauge and computational output set O CO mp- Then, the relations 
q = Ts + Hg mod 1, and o = Zs + Rg mod 2 for an optimal output o are unique. 

That is, once the resource state |\E'), the measurement planes and gauge? O comp are fixed, there 
is no freedom left to choose the classical processing relations. They are uniquely determined by the 
former. In particular, for fixed stabilizer state |^) and measurement planes, T = r(J gaU g e , O comp ), 

H — i?( /gauge; Ocomp) Gtc. 

A corollary of Theorem [l] is that given the measurement planes and an extremal pair /gauge , 
Ocomp; the resource state |^) uniquely determines the influence matrix T. One may ask how 
restrictive a condition the extremality of the pair /gauge, Ocomp is. In this regard, note 
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Theorem 2. Consider an MBQC on a fixed resource stabilizer state for fixed measurement planes, 
with an influence matrix T and input and output sets I(T), 0(T), such that no qubit a £ I c can 
be individually gauged with respect to I(T), 0(T). Then, there exists an extremal pair I gauge Q I, 
Ocomp O such that T = T(Ig aU g e , C omp)- 

The input set I(T) and the output set 0(T) which appear in Theorem [2] are uniquely specified 
by T through Definition [5j T(/ gauge , O com p) is uniquely specified by the pair /gauge, O com p through 
Theorem [TJ 

Theorem [2] states that all temporal relations for an MBQC, subject to the extra condition on 
the qubits which can be individually gauged, arise from extremal pairs /gauge, Ocomp- 

By establishing Theorem [2] we trade the condition of the pairs /gauge, O com p being extremal 
for the condition that no qubit in I c can be individually gauged. The latter is a more meaningful 
condition. Suppose a qubit a in I c could be individually gauged wrt /, O. Then K{a) exists. 
For any b with K(b)\ a = K{b) := K(b)K(a) is a valid correction operator for qubit 6, and 

K(b)\ a = l( a \ Hence, a could be removed from all forward cones and thereby be made a qubit 
in I. By imposing the extra condition in Theorem [2j we exclude temporal relations where certain 
qubits could be in the input set I but aren't. 

Theorem [l] is mute on the question of which extremal pairs /gauge, O com p are admissible. Theo- 
rem [3] below describes how much freedom remains for the choice of the classical processing relations, 
given C7(|¥)). 

Theorem 3. For MBQC with a fixed resource stabilizer state \^f) and fixed measurement planes, 
the classical processing relations for extremal I gauge , O comp , as specified by the matrices H, R, T , Z 
and the sets Igauge, O comp , o,re in one-to-one correspondence with the bases of the matroid Q(\fy)). 

After we have justified our restriction to extremal pairs of gauge input and computational 
output sets in Theorem [2] and have characterized the set of temporal relations compatible with a 
given resource state and set of measurement planes in Theorem [3j we now return to Theorem [TJ 
and show that a converse also holds. 

Theorem 4. Consider an MBQC on a stabilizer state with classical processing relations 

q = Ts+Hg mod 2, o = Zs+Rg mod 2 for an optimal classical output o, such that rkH = rkZ . 
Then the classical processing relations uniquely specify the stabilizer generator matrix Q{\ty)) in the 
a^/ag-basis, i.e. the resource stabilizer state |^) and set E of measurement planes up to equivalence. 

5.2 The MBQC - spacetime correspondence 

Our proposal of MBQC as a toy model for spacetime is inspired by Malament's theorem [40J, 
stating that for a continuous spacetime manifold, the metric is determined by the temporal order of 
spacetime events up to a conformal factor. Our Theorem [4] shows close resemblance with Malament's 
theorem if we assert the following correspondence: 



MBQC 



spacetime 



qubit 

measurement angles, 
one per qubit 

resource state | H 1 ) and 
measurement planes^ 



spacetime point 

conformal scale factors, 
one per spacetime point 

spacetime topology 

(metric up to conformal factor) 



(30) 
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While the first line of the correspondence (30) appears intuitive, the second line really comes about 
by identifying the left-over pieces on either side. Nevertheless, it seems a match: A reel-valued 
parameter, the measurement angle at a given qubit is mapped to a real-valued parameter, the scale 
factor at the corrsponding spacetime point. (Note, though, that the domain of the measurement 
angle is a compact set while that of the scale factor is not). A physical interpretation for this 
correspondence remains to be found. 



The support for the conjectured correspondence (30) comes from its third line. By analogy with 
Malament's theorem, the temporal order of measurement events should, up to equivalence at least, 
determine the resource state |^) and the set S of measurement planes. Theorem [4] almost gives 
this result. But the correspondence seems imperfect: It is not the temporal order alone — contained 
in the influence matrix T — which yields the pair |\&), £ up to equivalence. Rather it is the entire 



classical processing relations Eq. (15a), (15b), specified by the four matrices T, H, R, Z. 

But we can do better! Namely, it can be shown that the four matrices T, H, R, Z can 
be assembled to an influence matrix T cxt of an MBQC on a slightly bigger resource state \^'). 
Specifically, the support of \*$>') is, in comparison to \^f), enlarged at the temporal boundaries 
-^gauge> O CO mp- By this extension, all information contained in the classical processing relations 
becomes temporal information. 

To obtain T cxt we first choose a pair /gauge C I, O comp C O (1,0 follow from T), and extract 



from H,R,T,Z the corresponding matrices H, R, T, Z wrt to the normal form Eq. (26) specified by 
gauge , O comp . We assemble the composite matrix 



which is a square matrix of size (|0| + |/ ga uj 











°\ 


I 











H 


T 








\r 


Z 


I 


o / 


el + 


O c omp 


)x( 



gauge 



(31) 



| Ocomp | ) • 



T ext is the influence matrix for MBQC on a bigger resource state |\&') constructed from |^). 
For the support Q' of \^') we require two additional sets of qubits, I' and O' , with \I'\ = \0'\ = 

QUl'UO' . \$?') is obtained from \ty) by the following construction: 



'gauge | 



I O comp | , such that O' 




(32) 



Therein, the gates O — O = A s are cr s -controlled tr s -gates, i.e., Ascr^At = Og , A s er^A| = er^ (8) 
, etc, and the extra qubits in / and O are initially prepared in the eigenstate of as with 



eigenvalue 1. With the definition Eq. (32) of 1^'), the labelling of the blocks of rows and columns for 
the matrix on the r.h.s. of Eq. (65) is / / |(0 C omp) c |O comp |0 / for the columns and / / |/gauge|(-^gauge) c |0 / 
for the rows. 

We thus find that all information in the classical MBQC processing relations is temporal infor- 
mation for the computation on a slightly extended resource state. 
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5.3 Proofs of Theorems CDS] 

Theorem [T] is an immediate consequence of Lemma [3j 

Proof of Theorem^ Assume that I is valid input set and O is a valid output set for a given MBQC. 
Then, the stabilizer generator matrix of the resource state can be written in the <70/cx s -basis as 



I 
\ 



I 


r 


o c 








f T 


I 


A 


B 


c 





D 



(33) 



/ 



for some matrices A, B, C and D. The influence matrix T can be obtained as 






f T 









(34) 



with the column split between I and 1°, and the row split between O c and O. The matrices B 
and (B\C) do not necessarily have maximal row-rank. By row transformations of (-B|C||0|-D) we 
extract the dependent rows, and obtain 



V 



I 


r 


o c 








T T 


I 


A 











D' 





C" 





D" 


B'" 


C" 





D'" 



(35) 



/ 



However, note that each row in the third set of rows in the above matrix can either be interpreted 
as a correction operator K(a), with non-empty forward cone /c(a), for some a E O, or a gauge 
operator K(i) for some i G I c . The former is ruled out because every a £ O must have an empty 
forward cone. The latter is ruled out by the assumption that no qubit outside the set / can be 
individually gauged. Therefore that set of rows must identically vanish. 



Since ^(|^)) has full row rank, so does the matrix D' appearing in Eq. (35). We may then 



choose a set AO C O such that the columns of D' indexed by AO form a maximal independent set. 



We set O 



comp • " 



0\AO. Then, by further row transformations which do not affect T, the matrix 



in Eq. (35) can be converted to 



O c 



AO O, 



comp 



V 






fT 


I 





A' 











I 


A" 


B"' 


C" 








D"' 



(36) 



B'" has full row rank by construction. We can therefore find a set /gauge Q I such that the columns 
of B'" indexed by gauge form a maximal independent set. For any such set gauge we can convert 



the matrix in Eq. (36) fully into the normal form Eq. (26) without affecting T 
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For any of the above choices for /gauge Q I and O com p Q O, the resulting influence matrix 
/"(/gauge, Ocomp) can be extracted as 







T T 


/"(/gauge > O c omp) — 










Vo 






(37) 



with the column split between I and I c , and the row split between O c , AO and O CO mp = 0\AO. 
By comparison of Eqs. (34) and (37) we verify T = T(I glluge , O comp ). □ 



Remark: Comparing Eq.(36) with the normal form in Eq. (26), we can write the stabilizer 



matrix in a slightly varied form that can be useful later. 







/gauge 


AI 


r 


o c 


AO 


Ocomp 




/ 












T \ 












T T 


I 





















I 


Z T 






I 




Hi 








R T 














J 



(38) 



where Z = (Zi | Z2) and where E T = (Ej \ H^). 

Proof of Theorem^ Denote by £> the set of bases of Q, and by T the set of extremal classical 
processing relations of form Eq. (15a), (15b), specified by the triple (I gaU g e , Ocomp, {H, R, T, Z}). 
Then, the mapping h : B — > T exists and is a bijection. (1) Existence of h: By the normal form 
Eq. (26) of Q = (<^||5), for a given basis B(Q) the sets /gauge and O CO mp are extracted as follows. 
A qubit i is in /gauge if and only if the corresponding column of appears in B{Q). A qubit a is in 
Ocomp if and only if the corresponding column of S does not appear in B{Q). Knowing /gauge and 
Ocomp, {H, R, T, Z} is uniquely determined via Theorem [l] (2) Surjectivity of h: By definition of 



"extremal". (3) Injectivity of h: given /gauge and C omp> 



B{G) 



s o cornp ) is unique. □ 



Proof of Theorem^ We divide the proof into the following steps: i) First, with rk// = rkZ, the 
processing relations stem from an extremal pair /g aU ge;0 C omp- We show that given an extremal 
pair /gauge, Ocomp, the matrices H, R, T and Z are uniquely determined by the classical processing 



relations Eq. ( 15a|15b ). ii) From these matrices we can derive the corresponding normal form of 
the resource state uniquely. Let us denote this by M. In iii) and iv) we show that the following 
diagram commutes, which establishes the equivalence of normal forms of all extremal pairs. 



(/gauge; Ocomp) 
T,H,Z,R 

M 



M(A) 



-> (/' 



O' ) 

gauge "comp/ 



T',H',Z',R' 



(39) 



In Eq. (39), M(A) is a transformation on the normal form J\f, dependent on A. Because of their 
independence, we consider the transformations /gauge 



iii) and iv) respectively. 



^gauge and Ocomp 



O CO mp separately in 
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Extracting H, R, T and Z: Given any set of inputs I and outputs by Theorem [2j we know that 
there always exists an extremal pair (/gauge, O CO mp), where /gauge C /, O CO mp C O. By Eq. (29), 
T is uniquely specified by T. By the assumption of the classical output being optimal, Z can 



be brought into the unique normal form of Eq. ( 23 ) by left multiplication with an invertible 



matrix and (permuting the columns if necessary), Z can be extracted. The matrix H may not 

appear in Eq. (15a) in its normal form Eq. (27), nonetheless for an invertible matrix A, the 

-1, 



vector q specifying the measurement bases is invariant under H 



HA, 



A~ L g. By 



definition of / 
O 



gauge. 



every qubit in /gauge can be individually gauged with respect to / 



gauge. 



comp' 



Therefore, (up to row permutations), we can choose A such that 



HA 



(40) 



where the row split is between /gauge (upper) and (/ gaU ge) c (lower). A and H are unique. The 
classical output o is invariant under the transformation R — > RA' , g — > (A') _1 g. However, 
since Eqs. (15a) and (15b) refer to g in the same basis, A' is now fixed: A' = A. Then, R = RA, 



with the A of Eq. (40) 



Assembling the normal form: Using the unique matrices H, R, T, Z, by Lemma |3j we can now 
assemble the normal form Eq. (26) of the stabilizer generator matrix for the resource state 
|^). By assumption of Theorem 4j the matrices H, R, T, Z describe a valid computation, and 
the normal form Eq. ( 26 ) derived from them must thus yield a valid description of a quantum 



state. In particular, all Pauli operators specified by the rows of G({K, R, T, Z}) in Eq. (26) must 
commute. (The rows are independent by design of the normal form.) We have thus constructed 
a description of |^). Since H, R, T, Z are unique, so is \^f). 



We now proceed to construct the stabilizer 5(|\E')) from the processing relations when / f 



auge 



and Ocomp are not specified. From the classical processing relation Eq. ( 15a) we can still extract 



the input set / and the output set O, by testing which rows and columns of T identically vanish. 
Then, the possible choices for /gauge and O comp are limited to /gauge C / and O CO mp Q O by 
definition. 

Equivalence under A 

version of the normal form of G(\^f)) as shown in Eq. (38) which also includes additional 



'gauge 



/' 



gauge- 



In order to prove this we rely on the slightly variant 



detail about the correction operators. As noted above, the different normal forms for H can 
be interconverted by right-multiplication with an invertible matrix A (change of basis for g) , 
i.e. H, R — > HA,RA, where H T = (/|H T ). Under such a transformation, the upper part of 



the normal form Eq. (26) for Q{\^>)) remains unchanged, and the lower part is transformed 
(H T \ \0\R T ) — > A T (H 1 | |0|R T ). Invertible row transformations on Q(\^)) leave |^) unchanged, 
and 1^) is thus independent on the precise choice of /gauge C /• 



does not 



Equivalence under A : O com p —> 0' comp . Proving that a different choice of O CO mp 
change the stabilizer state is a little more complicated. We proceed in the following manner. 
From Theorem [3j we know that /gauge U Ocomp are the bases of a matroid. Therefore for 
two distinct computational output sets O com p and O c ' omp , there exists another computational 



output set 0' cc 



O c 



Therefore, 



',■>!(:]> ■ ^comp\{*}U{j}, where i £ Ocomp\O com p and j G O comp \O com p. 
it suffices if we show that the stabilizer state does not change if we change the computational 
output set from O com p to O' comp . Assume that the classical relation for the computational 
output set Ocomp is given as 



o = Zs + Rg = (Zi |Z 2 |/)s + Rg. 



(41) 



where the column split of Z is between 0c Omp , AO, and O CO mp- 



The corresponding normal 
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form (with the column split in cr^-part 7 S 



gauge | 1 gauge 



IS 






f T 


I 
















I 






E 1 ' 








R r 



(42) 



where T T 







Suppose that we transform O comp to O com p \ {i} U {i}, where i 6 O comp 



and j G AO. Without loss of generality assume that i is the last column of Z 2 . (It cannot be an 



all zero column because, then it would not be possible for it to be in O, 



comp' 



Let Z 



Z A 



and Z2 



O' 







Z B 


b 



. Then A 








b 





acting on Z achieves the transformation O com p to 



comp - 



A(Zi|Z 2 |/) 



x T 


T 

a 


1 


1 





Z A + bx T 


Z B + ba T 





b 







T 

a 


1 


1 





Z A + bx T 


Z B + ba r 


b 








(z'lll^P) 



(43) 
(44) 



We claim that this same transformation can be effected by row transformations of Q(\^)). 
First let us focus on the middle set of rows in Q(\^f)), namely the correction operators for AO. 



Then acting by M(A) 





a 





1 



gives us 



M(A)(0||0|I|zJ 









I 


a 





Zl + ab T 











1 


1 


b T 








/ 





a 


Z T B + ab T 











1 


1 


b T 



|o|/|zf). 



(45) 
(46) 



Now if take the last row in Eq. (46), namely ( 1 1 1 1 1 1 ) = c and add xc to the top set of rows 

(47) 



in Eq. (42) we obtain 

( I T T || I I I x I I Z T A + xb T ) ~ ( I f T || I I I Z/ T ) 



showing the equivalence of Z and Z'. The equivalence of R and R' under A can be shown in 

exactly the same fashion as for Z\ and Z^. 
This concludes the proof that the extremal classical relations completely determine the stabilizer 
state. □ 

6 Invar iance under the gauge transformations 
6.1 Gauge transformations and temporal order 

In this section we provide a different angle at Theorem [TJ namely we show that the gauge trans- 
formations Eq. ( 18 ) impose severe constraints on the possible temporal orders for a given stabilizer 



state and set of measurement planes. 



The gauge transformations Eq. ( 18 ) act on q and s, and can therefore have a non-trivial effect 



on the classical processing relation Eq. (15a), q = Ts + Hg. We now study this effect. Since the 



transformations Eq. (18) are caused by the insertion of stabilizer operators into the state overlap 
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in Eq. ([5]), they do not change the physical situation. Therefore, the temporal relations before 
and after any such transformation must be equally valid, although not necessarily identical. By 
insertion of the stabilizer operator into the overlap ($i oc |\l/), the stabilizer of \^>) and the the sets 
gauge, Ocomp do not change. Therefore, by Theorem [TJ the matrices T and H do not change. Thus, 
besides q and s, all that can change in the relation Eq. (15a) under a transformation Eq. (18) is g. 



The following two viewpoints are always equivalent: (A) The relation q = / g (s), under the action 



Eq. ( 18 ) of a Gk on (s, q) is changed into an equivalent such relation q = / g / (s) , with Gk '■ g — > g'- 
(B) The relation q = /(s,g) remains invariant under all transformations Gk, acting on the triple 
(s, q, g). We choose the latter viewpoint. 

We now infer the action of the transformations Gk on g. Without loss of generality we assume 



that the relations Eq. (15a) are given with H in its normal form Eq. (27), 



Ts 



mod 2. 



Furthermore, we assume that the pair /gauge) Ocomp is extremal. Since /gauge Q I by definition, 
for each i G /gauge, Qi only depends on g but not on the measurement outcomes s, q% = Qi- Now, 



the correction operators K(a), a G (O com p) c derived from the normal form Eq. (26) of Q(\^)), 
K( a ) I /gauge has no o^-part. Therefore, the corresponding transformations Gk do not flip qi, for all 
In order to preserve the relation qi = gt, they thus leave g unchanged. Now consider 



i G I, 



gauge- 



the other stabilizer generators, K(i), i G /gauge, obeying the conditions Eq. (22). By construction 



Gftu\ flips qi but no other qj, for i / j G /gauge- Hence, to preserve the relations qi = g%, 



flip 9ii but no other gj, i ^ j G /gauge- Thus, for a stabilizer element K — <S> a en 

G K ■ g — >g®w|r e . 



it must also 

(ah 



(48) 



Therein, we have assumed that the basis choice for g is such that the matrix H appearing in 



Eq. (15a) is of normal form Eq. (27). 



We have now fully specified the action of Gk on the triple (q, s,g), c.f. Eq. (18), (48). MBQCs 
satisfy the invariance condition 

q = Ts + J ?/g mod2^G K (<l)=TGK(s) + HGK(g) mod 2, \fK G «S(|tf». (49) 



It is evident that the requirement (49) of invariance of the processing relations (15a) under the 



gauge transformations poses constraints on the possible matrices T and H. In fact, as we show 
below, given C omp the matrices T and H are uniquely specified uniquely by the above invariance 
condition. 



To check the invariance condition Eq. (49) in a specific case, we return to our 3-qubit cluster 
state example of Section 2.2 We consider the effect of the transformations induced by generators 
Kl = o^o?\K 2 — 
earlier, / 



a^afaf^ and K 3 



(2) (3) 

<7s a , on the processing relations Eq. (16). As noted 



gauge 



{1}. Then, with Eqs. (18) and (48), 



Gk x 

Gk 2 
Gk 3 



q 
q 
q 



q® (1,0, 0) T , 
q® (0,1, 0) T , 
q® (0,0, 1) T , 



s© (0,1, 0) T , 
s©(l,0,l) T , 
s© (0,1, 0) T , 



,9i 
91 
9i 



9i 

91, 

91- 



1, 



(50) 



It is easily checked that the relation Eq. ( 16 ) is invariant under the transformations Gk 1 , Gk 2 an d 



Gk 3 of Eq. (50). However, if the transformations are restricted to q, s, the relation Eq. (16) is no 
longer invariant under the transformation induced by K±. 



We now return to the general case and show that, given the set Ocomp and the action Eq. (18) 



(48) of the gauge transformations on the triple (q, s,g), the invariance condition Eq. (49) uniquely 



specifies the classical processing relations Eq. (15a) for the adaption of measurement bases 
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Recall that we write the stabilizer generator matrix for |^) in the a^,/a s -basis as = 
($>\\S). Then, for the stabilizer generator K a £ S(\^f)) corresponding to the a-th row of Q(\^/)), 



with Eq. (18) the action of the gauge transformation Gx a on s, q is 

G Ka : s — > s row a (S), q — > q row a ($). 



With Eq. (48), the action of Gx a on g is 

G Ka :— > g©row a ($)k 



(51) 
(52) 



Here, row a ($)|/ gaugc denotes row a ($) restricted to /gauge- Then, the condition Eq. (49) for invariance 
of q = Ts + Hg under Gx a becomes 

row a ($) = Trow a (S) + #row a ($)| We mod 2. 

This condition must hold for all stabilizer generators K a simultaneously, hence 

mod 2. 



$ T = TS T + H<S> T \ L 



gauge X fl 



(53) 



By definition, the qubits in /gauge have empty backward cones, and the qubits in C omp have empty 
forward cones, hence T is of the form 








T 






where the column split is (0 C omp) c |O comp and the row split is /g aU ge|(/gauge) c , c.f. Eq. (29). By 
right-multiplication of relation Eq. (53) with a suitable matrix, we transform & T \i au cxn into a 
matrix of form (I|0) where the column split is between /gauge and /g aU ge- By definition of /gauge, 
such a transformation is always possible. Under the same transformation, 








$1 


$ 2 



s 



I (Ocomp ) c X 



(Si|5 2 ). 



(54) 



Inserting the above into Eq. (53), we find that H must be of normal form Eq. (27), H 
and 

$1 = TSi + H mod 2, 
$2 = TS2 mod 2. 



(55) 



Now, S2 must be an invertible matrix. This is the condition that, by definition of O CO mpj every 
measurement outcome in (O. 



comp J 



is correctable. Then, 



T = ^S^" 1 mod 2, H = $1 + <f> 2 S 2 ~ 1 Si mod 2. 
Hence the relation q = Ts + Hg is uniquely specified. 

6.2 Gauge transformations and computational output 



(56) 



In addition to Eq. 49 , we also require invariance of the classical output under the transformations 
o = Zs + Rg mod 2 = ZG K {s) + RG K {g) mod 2, \/K £ S(\^}). (57) 



Eq. (18), (48) 



Like Eq. (49), Eq. (57) is a determinism constraint. If for a single output bit o the relation 
o = z T s + r 1 g is not invariant under all gauge transformations Eq. (18), (48), then the value of o 



is guaranteed to be random, and thus useless as readout bit of a computation. Specifically, 
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Lemma 4. Assume an MBQC where the relation q = Ts + Hg is invariant under the gauge 



transformations Eq. (18), (48), but an output bit o exists whose defining relation o = z T s+ r 1 g 
is not invariant under the action of Gk for some K G 5(1^)). Then, the value of o is completely 
random, independent of the choice of measurement angles. 

Proof of Lemma^ For simplicity, consider first the special case where K £ 5(1^)) acts trivially 
on g, Gk{s) = g, V<?. We may then write o = X^iej s i + c f° r an offset c = r T g. We call the string 
s\j of measurement outcomes on J even (odd) if it has even (odd) weight. We denote the local 
post-measurement state on qubit a by \(p a , s a , q a (s, g)), where ip a is the measurement angle, s a the 
measurement outcome and q a specifies the chosen measurement basis. 

Under the transformation Gk, s — > s © As^, where, by assumption, As^| j is odd. Now, the 
probability of outputting o = c is 



p{o = c) 



E 

s| j=even 

E 

s| j=even 

E 

s| j=even 

E 

s| j=even 

E 

s| j=odd 
p(o = c). 



a{(Pa,S a ,qa(s,g)\ \$>) 



\aefl 



a{(Pa,S a ,qa(s,g)\ K\V) 



\aefl 



ci) qa 

®Aq K „\ |*) 



\aefl 



(58) 



ai^ai s a © Ask, ci) qa 

(sffi As K , S )\ |*) 



\aefl 



a{fa,Sa,q a (s,g)\ \V) 



Thus, p(o = c) = p(o = c) = 1/2. Note that in transitioning from the third to the fourth line of 



Eq. (58) we have used the invariance property Eq. (49), i.e., the assumption that the adaption of 



measurement bases is deterministic. 

In the general case, Gk '■ s — > s© Ask, g — > gffi Ag K - We note that we can choose any gauge 

2 



fixing g, and thus p(o = c) = | ^ ^ ( (^) a (<Ai, s a , q a (s, g)| j |* 

g s|j=cvcn \aef2 / 

analogous to the above we then find p(o = 0) = p{o = 1) = 1/2. □ 



By an argument 



7 What is time in MBQC? 

It has been established earlier that 'time' in measurement based quantum computation is a transi- 
tive binary relation. If a given MBQC is to be deterministically runnable, then we require a strict 
partial ordering among the measurement events. If, on the other side, in pursuing an analogy with 
general relativity we allow for closed time-like curves, then antisymmetry and irreflexivity are not 
required. Given this setting, in the previous sections of this paper we have analyzed the constraints 
on MBQC temporal relations imposed by a group of gauge transformations. 

In this section, we provide another angle at the question of what time in MBQC is. Namely, 
we introduce a second group of symmetry transformations on MBQCs which is related to local 
complementation [27J - |29j . These transformations in general change the MBQC on which they 
act, but leave temporal relations invariant. The possible temporal relations in MBQC thus label 
representations of this group. 



23 



7.1 Flipping measurement planes 

According to Eq. ([3]), for any qubit a in a given resource state, the local observable measured to 
drive the computation is 

O a [q a ] = cos (f a oT" + 



sm y? a cr 



Therein, q a is a linear function of the measurement outcomes {sb, b S 0}, c.f. Eq. (15a). 



Let's see what happens if we use a different rule for the adjustment of measurement bases, 
namely 

O'M = (-l)«»cosp a< #> +smip a a<$, (59) 

That is, if q a = 1, to obtain O r a [l] we are flipping the observable O a [0] about the <r s 0-axis rather 
than the 



-axis. Comparing Eqs. R3h and (59), we find that 



O'Ma] = (-l) ga O a [q a ] 



(60) 



independent of the measurement angle ip a - The measurements of 0' a [q a ] and O a [qa] are always in 
the same basis for the same q a , and the measured eigenvalues differ by a factor of (— l) qa . 

We call the transformation r s [i] : Ojjfft] — > O'^qi] flipping of the measurement plane at qubit i. 
On the elementary degrees of freedom, namely the resource state \^f), the Pauli observables a, 
a< sl (action on °s is implied), and the measurement angle ipi, the flipping r s [i] acts as 



S ) 



t, i : 



->■ cr 



-l)*f 



(61) 



The action of T s [i] on the Pauli operators a^' and the measurement angles ipj, for j ^ i is trivial. 
All our considerations are independent of the values of the measurement angles. In particular, the 



second part of the transformation Eq. (61) does not affect temporal order. 



We now discuss the effect of the flipping T s [a], a G Slona given MBQC. Denote the measurement 



outcome of a measurement of 0' a (q a ) by s' a . By Eq. (60), the following two measurement procedures 
are always equivalent. (I) Measuring O a [(/ a ] and outputting s a , and (II) Measuring 0' a [q a ] and 
outputting s' a + q a mod 2. We may call the device that performs Procedure I a 0-box, and the 
device that performs Procedure II a s<ft'-box. Then, 





(|)-box 



s4>'-box 



The prime in the s^'-box accounts for the fact that not the measurement outcome s' a itself is 
outputted, but rather the locally post-processed value s' a + q a mod 2. Now, instead of outputting 
s' a + q a , the device at a may only output s' a (that is an s^-box), and the classical post-processing 
relations for the adaption of measurement bases are modified accordingly, i.e. s a — > s' a +q a mod 2. 
Can the resulting relations again be written in a form q = TV + H'g? 

We now attempt transforming a (p-box into an s(p-box at qubit a. The vectors of measurement 



outcomes s and s' are related via s = s' + e a e^q. Inserting this relation into Eq. (15a), we obtain 



(/ + Te a e^)q = Ts' + Hg mod 2. 



(62) 
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Case I: T aa = 0. Physically, this means that the measurement basis at the flipped qubit a does 
not depend on the measurement outcome at a, before the transformation. Multiplying Eq. (6i 
with from the left yields q a = e^Ts' + e^Hg mod 2. Inserting back into Eq. (62), we obtain 

T' = T + Te„eTT mod 2. 



Likewise, 



H' = H® Te a eiH, Z' = Z © Ze a e T a T, R' = R 



Ze a e T a H. 



(63) 



(64) 



Eqs. (63) and (64) completely describe the effect of the flipping T s [a] of the measurement plane at 



qubit a on the classical processing relations Eq. (15a), (15b). 



Remark: If the matrices H, R, T, Z are given their normal form Eq. ( 27 ) wrt the pair /gauge , O com p 
then the flipping of the measurement plane at any vertex a with T aa = leaves this normal form 



intact. We can therefore state a transformation rule equivalent to Eqs. (63), (64) for the matrices 
H, R, T, Z. This rule is, in fact, simpler. We again consider the composite matrix T ext , c.f. Eq. (31), 



( 










I 











H 


T 








u 


Z 


I 


o / 



The effect of flipping of the measurement plane at a then is 

r ext — > T e ' xt = T cxt + T C xt e a ej T e x t mod 2. 



(65) 



Thus, the rule is just the same as Eq. (63) for the original influence matrix T. 

Case 2: T aa = 1. In this case, the measurement basis at a does depend on the measurement 
outcome at a. This is an example for a closed time- like curve (only involving the measurement 
device at a), and an obstacle to deterministic runnability. Now, the matrix / © Te a eJ on the left 
side in Eq. (62) is not invertible, e^(I © Te a e^) = 0. Hence, the relation (62) can not be solved 
for q in this case. There is no relation Eq. (15a) with the same sets -/gauge; Ocomp before and after 
flipping. 

We now discuss the consequences of flipping measurement planes for the above two cases. 
7.2 Flipping measurement planes and local complementation 

We now return to the above Case 1, namely when flipping of a measurement plane yields a compu- 



tation with a new relation Eq. (15a). Note that the computation before and after the flip generate 
the same output distribution. Flipping a </>-box into an s^'-box is an equivalence transformation, 



only based on the operator identity Eq. (60). Changing an s^'-box into an s^i-box is again an 
equivalence transformation, provided it can be carried out. 

The influence matrices T and T' before and after the flipping, respectively, are in general not 



the same, c.f. Eq. (63). However, T and T 1 still generate the same temporal order, as we now show. 

Lemma 5. Be T an influence matrix with T„ = 0. Then, T and T' = T © Te^e^T generate the 
same temporal relation under transitivity. 



Proof of Lemma [5| Let's introduce a shorthand a — > c for c G fc(a) (meaning that the 
measurement outcome at a influences the measurement basis at c). Now, we have to show that 
e <t f < ^=^ > e G^t' f, for any T' generated from T by the transformation Eq. (63). 
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(I) "= 

mi — > m,2 



Assume that e -<t f ■ Then there exists a sequence of measurement events e 



rewrite the transformation rule Eq. (63) for the flip T s [i] as 



/. Can we break the arrow a 
i] 



c, say? To investigate this, let us 



t s \i 



fc(a) 
fc(a) 



fc(a) € 
fc(a), 



fed) 



if i £ fc(a) 
if i g fc(a) 



(66) 



Case 1: a — > i before the transformation T s [i]. Then, fc'(a) = fc(a) © fc(i). Since Ta = by 
assumption, a — > i after the transformation T s [i]. Sub-case la: i c before the transformation 
T s [i\. Since Ta = {i fc(i)), i — > c after the transformation r s [i]. Thus a — >■ i — >• c after the 
transformation, and hence a -<y/ c. Sub-case lb: i -/> c before r a [i]. Then, a — > c remains after the 
transformation. Case 2: a i before the transformation t s \i\. Then a — > c after r s [i]. Thus, in all 
cases a ~<x' c -> an( i therefore e ~<t' f ■ 
(II) " < j=" : From Eq. & r,[i] 



Apply a series of transformations Eq. (63) on an initial influence matrix T with vanishing 
diagonal part may produce an influence matrix with a non- vanishing diagonal part. Thus, the 



application of the transformation Eq. (63) is restricted. To circumvent this problem, we introduce 
a modified transformation 



f[»] : T — ► T = T + Te;ef T + P(Te 4 ef T) mod 2. 



(67) 



Clearly, this transformation takes influence matrices with vanishing diagonal part to influence ma- 
trices with vanishing diagonal part, and thereby avoids the problem of restricted applicability of 



transformation Eq. (63). Note that the transformation Eq. (67) has the form of local complemen- 



tation, albeit the influence matrix T that it acts on will in general not be symmetric. 



But what is the physical significance of transformation Eq. (67)? The only additional effect of 
the transformation f[i] 
the transformation, 

exchanges cr S( f, O c s on a respective qubit. The action of f[i] on the elementary degrees of freedom 
therefore is 



is the cancelling of the diagonal part of the influence matrix after 
This can be achieved by a local unitary that 



over t. 

c.f. the last term in Eq. (67) 



(0 



Tl 







<pi - 


(-1) 




sq> ' 


1*) - 





ri 



Vj € fcH) n 6c(»), 



(68) 



We find that the local measured operators for all qubits j G 6c(i) D /c(i) change in a way that 
cannot be accommodated by a change of the respective measurement angle. For those qubits, 
the new measured observables lie in a different equatorial plane of the Bloch sphere. Therefore, 



the transformation Eq. ( 68 ) , unlike the transformation Eq. ( 63 ) , does not necessarily map a given 
computation onto itself. What it does, however, is mapping a given computation to a computation 
with the same temporal relation. 



Lemma 6. Be T an influence matrix with Ta = 0. Then, T and T' 
generate the same temporal relation under transitivity. 



V(Te ie fT) 



Proof of Lemma [6} Assuming the initial influence matrix has vanishing diagonal part, we split the 
transformation T — > T © Te^ef T © D(Tejef T) into two steps, namely T — > T = T © Te { efT 
and T — > T" = T © V(T'). By Lemma [5j T and T' generate the same temporal order. Now 
= 1 for some k € £1, k ^ i. This requires that Ty. = = 1. Then, we also 
1 and T" k = = 1. Thus, k -<t' k and k <t" k. The closed time-like curve 
involving k is not changed by setting T^ k = 0. All other relations trivially remain unaffected by 
the transformation T' — > T" . □ 



assume that T' kk 



have T[ k 



rpl 
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We may now consider the orbit Ot of a particular influence matrix T under local complementa- 
tion Eq. (68). Every local complementation f[i] permutes the elements of Ot- Therefore, we have 
a homomorphism f[i] — > P[i] where P[i] is a permutation matrix of size \ Ot\ x \Ot\- The matrices 
{•P[*]|i £ ^} generate a representation of the local complementation group. The representation 
acts on a given orbit Ot of influence matrices, and all T £ Ot give rise to the same temporal rela- 
tion under transitivity. Therefore, this temporal relation is a label of the given representation. All 
possible temporal relations in MBQC arise as labels of representations of the local complementation 
group. 



8 Breaking up closed time-like curves 
8.1 Closed time-like curves of length 1 

In this section we return to the case of Ta = 1 before flipping the measurement plane at qubit i. 
Such a qubit i cannot be in /gauge; since /gauge — I by definition. If Ta — 1 then i £ bc(i). The 
backward cones of all qubits in i are empty by definition of /, however. Likewise, i Ocomp- If 
Tjj = 1 then i £ fc(i). However, fc(a) = for all a £ O comp . Thus, there is only one case to 
consider, namely i G (/ g augc) c n (O comp ) c . 

In this case, there exists a correction operator K(i) for qubit i before the flipping, K (i) = 
a s< p <g> K{i)\q\s. After flipping at i, this operator turns into 

T 8 [t\(K(i)) = af ® K(i)\ a \i =: K'(i). (69) 



That is, the operator T s [i](K(i)) resulting from flipping at i is a gauge type operator, c.f. Eq. (22). 
Thus, the flipping transformation T s [i] (when Ta = 1) enlarges /gauge by one qubit, 

Ts[i] • /gauge ^ /gauge U if T%i — 1. 

Furthermore, after the flipping at i there no longer is a correction operation for qubit i, hence 

T s [i\ . Ocomp ^ Ocomp U if T{{ — 1. 

This has two consequences. First, the forward cone of i becomes empty. In particular Ta = after 
the flipping. Thus, the closed time-like curve consisting of qubit i has been removed. Second, an 
additional bit of optimal classical output is being created by the flipping at i. What does that 
output bit signify? 

Recall that before the flipping at i, the rule for adjusting the measurement basis at i is 



! 

= s 



i + Sj mod 2, 



for some set J C 17. Here, we have dropped a constant offset h T g on the r.h.s. The symbol "!" 
above the equality means that equality is a requirement for the correctness of the computation, but 
it cannot be deterministically implemented. As follows from Eq. (60), the measurement outcomes 
before and after the flip, Sj and s[ are related via S{ = s[ © qi. For all the other qubits, s'j = Sj. 
Substituting this into the above relation, we obtain 

s- + s'j mod 2 = 0, Vq' a £ Z 2 . (70) 

jeJ\i 

Thus, the additional output bit Oj = s' { + Y^jeJXi s 'j m °d 2 is a flag. If Oj = then the computation 
succeeded, and if o% = 1 then it did not. 
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Figure 3: Breaking a closed time- like curve of length 1. The looped qubit becomes an element of 
gauge H Ocomp after flipping the measurement plane. As such both its forward cone (0 C omp Q O), 
and backward cone (/gauge Q I) must be empty. 



Now suppose that the problem solved by the given MBQC is in NP. Then, this flag bit is not 
necessary. The remaining output may be efficiently checked for correctness anyway. Thus, one may 
safely discard the extra bit ot of output. Not post-selecting on o, = can, if anything, only increase 
the success probability of the computation. We thus arrive at 

Lemma 7. Be Mi an MBQC with a classical output o and influence matrix T such that Ta = 1, 
i.e., Mi has a closed time-like curve involving a single qubit i £ !1, Be Mi the MBQC with the 
same classical output o, obtained from Mi by flipping the measurement plane at i. Then, the 
closed time-like curve of i in Mi is removed in M2- Furthermore, if Mi solves a problem in the 
complexity class NP with probability p, then M2 solves the same problem with probability > p. 

Remark: Lemma [7] does not guard against the inefficiencies of post-selection, in particular if 
multiple CTCs of length 1 are being removed. While the success probability after removing the 
CTCs is guaranteed not to be smaller than for the original computation with the CTCs (which can 
only be executed using post-selection), neither it is provably significantly larger. 

Event horizons. Let us consider the flow of information between qubit i whose measurement 
plane has been flipped and the other qubits. Before the flip (MBQC Mi of Lemma i £ 
(O comp ) c n (/ g auge) c - After the the flip (MBQC M2 of LemmaJT]), i G C omp D -fgauge- In M2, since 
i S -fgauge, no information for the adaption of measurement basis is flowing into site i from the other 
sites. Likewise, since i £ O com pi no information for the adaption of measurement bases is flowing 



out of site i. Finally, because of the normal form Eq. (27), the measurement outcome appears 
in only one readout bit. This readout bit is Oj as given in l.h.s. of Eq. (70), which is precisely the 
bit of classical output that can be discarded if the problem solved by the quantum computation is 
in NP. If Oi is discarded, then no information is flowing out of the site i at all. Thus, in summary, 
from the viewpoint of classical processing, qubit i in M2 becomes entirely disconnected from the 
computation. It vanishes behind the MBQC counterpart of an event horizon. 

8.2 Closed time-like curves of length > 2 

A closed time-like curve C = {1,2, ..,/} of length I > 2 is a set of qubits such that the relations 
1 — > 2, 2 — > 3, .. ,i — > i + 1, .., I — 1 — > I, I — > 1 hold. We call a closed time-like curve C 
minimal if no proper subset of C is a closed time- like curve. Minimal closed time-like curves are of 
the form 

fc{i) n C = {i + 1}, l<i<l, (where 1 + 1 = 1). (71) 
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That is, for a minimal closed time-like curve, no arrows other than the defining ones are present. 
Proof: Assume 3k > 1 such that i + k G fc(i). Then, £ = {1, 2, .., i — 1, i, i + k, i + k + 1, .., 1} is 
also a closed time-like curve, and a proper subset of C Hence, C is not minimal. Contradiction. □ 
Regarding the embedding of C in the overall temporal structure, we have 



gauge; 



: n(o, 



comp ; 



c /. 



(72) 

Furthermore, Vi € £, 



Proof: Vi G £, i — 1 G bc(i), hence bc(i) ^ 0. Hence, i g" -fgauge 
i + le fc(i). Hence, fc(i) / 0. Hence, i g" O com p C O. □ 

Now, the CTC C can be broken up between qubits 1 and 1. This time, no transformation t s is 
necessary. The starting point "1" for the labelling is of course arbitrary. To break up the CTC, we 
enlarge the sets of gauge input and computational output, 



gauge 



^gaugc U {1}, O comp 



o, 



comp 



U{1}, 



(73) 



and modify the correction and gauge operators 
K(l) - 
K{a) - 

k{%) - 



K\l) 


:=K(l), 


K'{a) 


. J K(a), 
■ \ K(a)K(l) 


K\£) 


f m, 


\ K(i)K(l), 



if 1 fc(a) 
if 1 G, 

if K(i) 
if K(i) 



,Va G (O, 



comp, 



Z \{1}, 



(74) 



(i) 



,Vi G /, 



a 



gauge- 



The first line in Eq. (74) says that K(l) is re-interpreted as K{1). Hence, qubit I no longer has a 
correction operation and thus becomes a member of O CO mpj as required in Eq 
line in Eq. (74) ensures that K'(a)\\ = 1^ for all a G (0 C omp 
C 



i g r 

gauge 



(73) 



5 ) c such that 6c'(z) 



The second 
= 0. Thus, 



r 



as required in Eq. (73). The third line in Eq. (74) makes the gauge operators 
compatible with the new normal formfor {Jd 1 ! 7 ) (based on Ig augc and O comp ), and has no effect 
on temporal order. Since, after the reshuffling of the gauge and correction operators according to 



Eq. (74), 1 G -fgauge and I G O' comp and we have successfully broken up the CTC C between 1 and 



I. Since the computational output set has been enlarged by one qubit, /, we have one bit o\ of 
additional classical output. What does it signify? - When retracing the transformations Eq. (74) 



in the normal form Eq. (26) of the stabilizer generator matrix ^(|^)), we find that the matrix Z is 
updated according to 

bc(lf\{0 



Z\{1} 



(75) 



where bc(l)\{/} is the characteristic vector of bc(l), restricted to O C omp\{0> and Z\{1} is the 



matrix obtained from Z by deleting column I. Thus we find that the new output bit o\ is 

s a mod 2. 



o', 



81+ £ 

a£bc(l)\{l} 



(76) 



We compare this to the relation q\ = s; + Ylaebc(i)\{i} s « mod 2 before the transformation Eq. (74), 
and find the following interpretation of o\: The temporal relations before the transformation Eq. ( 74 ) 
contain a closed time-like curve and can therefore only be implemented probabilistically. So one 
may as well assume q\ = 0, perform the measurement and later check whether the relation between 
qi and s was obeyed. With Eq. (76), this check amounts to = 0, and d l is thus a flag for the 
correctness of the computation. If the computational problem solved is in NP then we can afford 
to discard this extra bit o[ of output. The solution can be efficiently checked anyway. In no longer 
post-selecting on o[ = 0, we retain all the 'good' cases while admitting additional ones. The success 
probability of the algorithm does not decrease. We thus arrive at 
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Lemma 8. Be M\ an MBQC with a classical output o and influence matrix T giving rise to a 
minimal closed time-like curve C = {1,2,..,/} of length > 2. Be M.2 the MBQC with the same 
classical output o, obtained from A4i by enlarging Igauge — > Igauge U {1}, O comp — > O comp U {I}. 
Then, the closed time-like curve C is not present in M.2- Furthermore, if Mi solves a problem in 
the complexity class NP with probability p, then .A/f 2 solves the same problem with probability > p. 



We still need to check that by the transformation Eq. ( 74 ) of the correction and gauge operators 



we do not create additional closed time-like curves in exchange for removing the minimal CTC C. 



In this regard, note that the middle line of Eq. (74) amounts to 



fc(a) — ► fc'(a) = fc(a) /c(Z),Va G (0' comp ) c with 1 G /c(a). 

Therefore, for any b G (^g aU ge) C > & e / c '( a ) om y if & G f c ( a ) or 6 G fc(l). Thus, b G fc'(a) only if 
a -< b before the enlargement of igauge, O comp . Therefore, no new closed time-like curves are being 
created. 



8.3 Removal of closed time-like curves 

Combining the above two cases of CTCs of length 1 and > 2, respectively, we can remove CTCs 
repeatedly until none remains. The process must stop, since at some point every qubit is in the 
gauge input or/and the computational output set. We thus arrive at the following 

Theorem 5. Be Mi an MBQC solving a problem in the complexity class NP with probability p. 
Then, all closed time-like curves can be removed from the influence matrix T of Mi, resulting in a 
deterministically runnable MBQC M2 which solves the same problem with probability > p. 

Remark: Theorem [5] does not guard against the inefficiencies of post-selection. While the 
success probability after removing the CTCs is guaranteed not to be smaller than for the original 
computation with the CTCs (which can only be executed using post-selection), we cannot prove 
that it is significantly larger. 



9 Conclusions and outlook 

In this paper we have suggested MBQC as a toy model for spacetime emerging from quantum me- 
chanics. We have analyzed the constraints on MBQC temporal order that arise from the quantum 
mechanical randomness of measurement, in combination with the computational requirement of pre- 
venting this randomness from affecting the quantum information processing. We have identified two 
groups of symmetry transformations, a gauge transformation which leaves every MBQC — including 
the temporal relations — invariant, and a second transformation, related to local complementation, 
which only preserves the temporal relations. We have shown that for any given resource stabilizer 
state \^) and set E of measurement planes, the compatible temporal relations all arise from the 
bases of a matroid encompassing I^^X. Finally, we have identified event horizons as a further piece 
of the phenomenology of general relativity which has a counterpart in MBQC. 
At this point, we are led to ask the following questions: 



1. We introduced a group of gauge transformations Eq. ( 18|48 ), and a group symmetry transfor- 



mations Eq. (67), generated by flipping measurement planes. Both transformations preserve 



MBQC temporal orders. Can the two groups be unified? 

2. Some of the temporal relations admitted by the matroid ^(|^)) contain closed time- like 
curves. Given a stabilizer state |^) and set of measurement planes S, can we find a similar 
algebraic (or other) structure which comprises only the partial orders? Can the partial order 
of measurements with the smallest set Ocomp be efficiently computed? 
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3. We have established that all temporal relations free of self-loops appearing in MBQC arise 
as representations of the local complementation group. However, these representations are in 
general reducible. If we decompose them into irreps, where does temporal order go? 

4. In proposing MBQC as a toy model for spacetime, here we have focussed on the temporal 
part. What is a suitable notion of space that can be associated with a resource state and a 
set of measurement planes in MBQC? 

5. In MBQC, the link between the randomness in quantum mechanical measurement and tem- 
poral order is the principle that the randomness of measurement outcomes should not affect 
the logical processing. In a context more general than quantum computation, what could this 
principle be replaced by? 

Remark regarding Question up With Theorem [4] we arrive at a situation closely resembling Mala- 
ment's theorem |40| . The latter states that for a continuous spacetime manifold, the temporal order 
of spacetime events determines the spacetime metric up to a conformal factor, i.e., determines the 
topology of spacetime. In MBQC, the classical processing relations — which are essentially tempo- 
ral information encoded in the influence matrix T, plus the extra matrices H, R, Z related to the 
boundary sets /gauge and O comp — determine the entire computation up to the measurement angles. 
To further pursue an analogy with Malament's theorem it is therefore desirable to have a charac- 
terization of 'space' in MBQC that admits a discussion in terms of topology. Cellular complexes 
thus seem a suitable choice. Indeed, certain stabilizer states considered as resource states admit the 
characterization in terms of a cellular complex, such as the planar code state |41j (two-dimensional 
complex) as well as computationally universal two-dimensional cluster states (three-dimensional 
and four-dimensional complexes). 

Remark regarding Question Uk A possible criticism of MBQC as a toy model for a quantum 



spacetime is that the laws which govern it, namely the classical processing relations Eq. (15a) 



(15b) are not laws of nature, but only rules imposed by the requirement of shielding the processed 
quantum information from the randomness of measurement. As such, the processing relations may 
be obeyed or violated at will by an operator running the MBQC. 

To this we respond that if the operating unit for a measurement-based quantum computer was 
a dedicated device, it would not have the freedom to violate the classical processing rules. For it, 
they would be the laws of nature. They would follow straight from Newton's axioms if the device 
was mechanical, and from Maxwell's equations if it was electrical. The processing rules could only 
be violated by a conscious being, such as a human or trained animal. 

Such beings, bound to the pull of gravity and the consequences of the no-cloning theorem by their 
belonging to this 'real' universe, but entitled by their free will to disobey man-made regulations, can 
violate the processing rules of MBQC precisely because they have an existence outside it. Therefore, 
a theory of MBQC should not be required to describe them. Neither quantum mechanics nor the 
theory of gravity make statements about hypothetical objects that jump in and out of spacetime. 

To conclude, we would like to recall the main idea underlying this work. In attempts to unify the 
theory of general relativity with quantum mechanics, one may take the viewpoint that spacetime is 
not an independent construct, but rather a consequence of the laws of quantum mechanics. Once 
this assertion is spelled out, the natural next step is to identify the key quantum property which 
yields a mechanism for generating temporal order, and to illustrate this mechanism in a toy model. 
Measurement based quantum computation provides such a toy model. Therein, the key property 
which leads to temporal order is the inherent randomness of quantum measurement. 
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